arXiv: 1501.02468vl [cond-mat.stat-mech] 11 Jan 2015 


Collective Motion with Anticipation: Flocking, Spinning, and Swarming 
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We investigate the collective dynamics of self-propelled particles able to probe and anticipate the 
orientation of their neighbors. We show that a simple anticipation strategy hinders the emergence 
of homogeneous flocking patterns. Yet, anticipation promotes two other forms of self-organization: 
collective spinning and swarming. In the spinning phase, all particles follow synchronous circular 
orbits, while in the swarming phase, the population condensates into a single compact swarm that 
cruises coherently without requiring any cohesive interactions. We quantitatively characterize and 
rationalize these phases of polar active matter and discuss potential applications to the design of 
swarming robots. 
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I. INTRODUCTION 

Over the last 20 years physicists have devoted signif¬ 
icant efforts to elucidate the numerous dynamical pat¬ 
terns observed in populations of living organisms. Flock¬ 
ing patterns are a prominent example. They refer to the 
self-organization of an ensemble of motile individuals into 
a homogeneous group undergoing coherent directed mo¬ 
tion. As initially suggested by Vicsek et al. [T], simple 
alignment interactions are sufficient to trigger the emer¬ 
gence of flocking patterns. Alinement-induced flocks have 
been demonstrated theoretically, numerically and in syn¬ 
thetic experiments m?]. In addition, quantitative anal¬ 
ysis of animal trajectories [SHTT] support that the large 
scale features of some animal flocks can be rationalized 
by simple behavioral rules at the individual level, includ¬ 
ing velocity-alignment interactions. However, inferring a 
model from actual data requires a minimal set of hypoth¬ 
esis on the possible form of the interactions m, or on the 
statistical properties of the observables [H]. Until now, 
the overwhelming majority of the available models has 
been merely restrained to couplings between the instan¬ 
taneous positions and orientations of neighboring individ¬ 
uals (see unumn] for noteworthy exceptions). Inves¬ 
tigating alternative dynamical rules may provide further 
insight into the collective dynamics of motile individuals. 

Here, we generalize the conventional description of po¬ 
lar active matter [5] . We consider the dynamics of motile 
particles, which aligns their velocity with the orientation 
of their neighbors by anticipating their rotational mo¬ 
tion. Naively, anticipation would be expected to yield 
more robust flocks. In striking contrast, we show that 
simple anticipation rules result in much richer collective 
behaviors: Firstly, when motile particles strongly antic¬ 
ipate the orientational changes, the population simulta¬ 
neously breaks a continuous and a discrete symmetry. 
The system self-organizes into a spinning state where all 
the particles follow closed circular orbits in a synchro¬ 
nized fashion. Secondly, a combination of alignment and 
moderate anticipation results in the emergence of a polar- 


liquid phase akin to the flocking pattern observed with¬ 
out anticipation. Finally, at the onset of the flocking-to- 
spinning transition, the population forms stable compact 
swarms despite the absence of any attractive couplings. 

The paper is organized as follows: we first introduce 
the equations of motion of the interacting self-propelled 
particles. The alignement and anticipation rules are de¬ 
scribed. We then discuss the emergence of the three or¬ 
dered phases (spinning, flocking, and swarming phases). 
They are quantitatively characterized, explained and 
compared to the conventional patterns of polar active 
matter. We close this paper from an engineering per¬ 
spective. We show that minimal anticipation rules can 
provide an effective and safe design strategy to interrupt 
the directed motion of a flock without shutting down the 
propulsion mechanism at the individual level. 

II. MODEL: ALIGNMENT AND 
ANTICIPATION 

We consider an ensemble of N self-propelled particles 
in a two-dimensional space. The particle i, located at 
r/t), moves at a speed vq along the unit vector pi that 
makes an angle 0{pi) = 9i with the a;-axis 

Yi{t)=VoPi. ( 1 ) 

The equation of motion of their instantaneous orientation 
defines their anticipation and alignement rules: 

Oi{t) = -^( sin [Oi - [Oj + acTj)] ) (2) 

where a is a scalar parameter and aj = 9j/\6j \ is the sign 
of the angular velocity. The particles have a finite inter¬ 
action radius R, and Uli denotes the ensemble of particles 
interacting with the particle. We henceforth refer to 
CTj as the spin of the particle j. This quantity is akin 
to the spin variable defined in |14j up to a multiplicative 
factor, which is the local curvature of the particle trajec¬ 
tory. In Eq. the angular noises ^/s are uncorrelated 
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Gaussian random variables of unit variance, and r is an 
orientational relaxation time. For the sake of simplicity, 
units are chosen so that r = 1 and i? = 1 , and the control 
parameters left are N, vq, a, r], and the system size L. 

Equations and have a simple physical meaning. 
When 0 = 0, they reduce to the continuous-time ver¬ 
sion of the seminal Vicsek model m- the particles in¬ 
teract via effective torques that promote alignment with 
the instantaneous orientation of the neighboring parti¬ 
cles. Note that neither the momentum nor the angular 
momentum is conserved in Eqs. 1 and 2. When a > 0, 
the mean orientation of the neighbors is anticipated in 
a simple fashion. If the orientation of the j**' particle 
rotates in the clockwise (resp. anti-clockwise) direction, 
an effective torque promotes the alignment along the di¬ 
rection 0j — a (resp. 0j -I- a), a is a constant angle used 
to anticipate the future orientation of the neighboring 
particles. 

In order to gain a better insight into this interaction 
scheme, we expand the sine functions and use elementary 
algebra to recast Eq. into 


0i{t) = — cos a ( sin ( 0 ^ — 0j) ) 




— sin a ( sin 0i — {0j + ^j~^) ) 


( 3 ) 

j(zn, + 


which is the linear superposition of two models: a 
continuous-time Vicsek model and a a = ^ model. The 
second term on the r.h.s of Eq. promotes alignment 
with the acceleration of the surrounding particles. As the 
magnitude of the translational velocity is constant, the 
acceleration of a particle is perpendicular to its direction 
of motion: 0{dtPi) = 0i + ai^. Equationj^is qualitatively 
similar to the discrete time model introduced in [13] . We 
revisit here the impact of the velocity-acceleration cou¬ 
pling from a different perspective and further character¬ 
ize the resulting phase behavior. 

The special status of the a = 0 and a = f models 
prompts us to characterize them separately. However, a 
comprehensive numerical characterization would require 
to vary systematically three independent parameters that 
are the propulsion speed vq, the particle density and the 
noise magnitude 77 . In this paper we set vq = 1/2, which 
is a conventional value that facilitates comparisons with 
previous polar active matter models [mill- If not men¬ 
tioned otherwise, we provide results corresponding to a 
particle density p = 8, corresponding to N = 2048 parti¬ 
cles. Smaller density values yield qualitatively identical 
results, but require longer convergence times for finite 
values of a. We therefore focus solely on the impact of 
the anticipation angle a and the angular noise rj. 

When a 7 ^ 0 Eq. 2 and 3 are implicit as di depends on 
the time derivative of 0i. We solve Eqs. andnumeri¬ 
cally using periodic boundary conditions in a square sim¬ 
ulation box of size L and a forward Euler method. This 
method requires specifying the discretization scheme for 
the spin variables. Aiming at describing the impact of 
anticipation, we have naturally defined the discretized 
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FIG. 1: Color online. Continuous-time Vicsek model (a = 0). 
(a) Superimposed plot of subsequent particle positions in the 
isotropic phase (77 = 0.5). The color accounts for the spin 
value. Red: cr = — 1, and black: cr = -|-1. (b) Instantaneous 
orientations and positions of the particles in the isotropic 
phase (77 = 0.5). To improve clarity, only one half of all 
particles is shown, (c) and (d) Same plots for the flocking 
state (77 = 10“^). (e) Polarization plotted as a function of 
the magnitude of the angular noise. For all the pictures and 
plots: p = 8 , 170 = 0.5 and L = 16 (A = 2048). 




model by setting (Ji{t) = sign[ 6 >i(t) — 0i{t — 5t)], where 
the time step 6t is 10 “^r in all the simulations. 


III. RESULTS AND DISCUSSION 
A. Case a — 0: Alignment-induced flocking. 

Here, we briefly recall the phenomenology of the model 
when a = 0 , which is equivalent to a continuous-time ver¬ 
sion of the classical Vicsek model m- For high noise val¬ 
ues, the self-propelled particles undergo uncorrelated per¬ 
sistent random walks. The population forms an isotropic 
and homogeneous gas depicted in Figs. b. Decreasing 
the noise value below 770 = 0.45, the rotational symme¬ 
try of the particle orientations is spontaneously broken: a 
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macroscopic fraction of the population propel in the same 
direction along straight trajectories, see Figs. [^,d, and 
Supplementary Movie 1) [TH]. This flocking transition, 
is quantitatively captured by the sharp variation of the 
order parameter If = |(pi)i,t| (the mean polarization), in 
Fig. [ 5 . 

Note however that, unlike what is found in simu¬ 
lations of discrete time Vicsek models (in the dilute 
limit) [min], we do not observe the propagation of band- 
shape excitations at the onset of collective motion. The 
reason for this discrepancy is purely technical. Based 
on previous numerical observations m using an angu¬ 
lar noise as in Eq. the emergence of localized bands is 
expected to occur only in extremely large systems. The 
relatively small size of our largest simulations {L = 16, 
N = 2048) explains why solitonic bands do not form here 
at the onset of collective motion. 

The simulations in the case a = 0 reproduce the salient 
features of the Vicsek flocking transition, and therefore 
validate our numerical scheme. 


B. Case a = Anticipation-induced spinning. 

We now investigate the anticipation model correspond¬ 
ing to a = Not surprisingly, in the high noise regime, 
all the particles behave like uncorrelated persistent ran¬ 
dom walkers again, and form a homogeneous isotropic 
phase with no preferred spin. Figs. [T}i and How¬ 
ever, reducing the noise amplitude below rjK = 7 x 10“^, 
two symmetries are spontaneously broken: the continu¬ 
ous rotational symmetry of the particle orientation, to¬ 
gether with the discrete Z 2 symmetry of the spin vari¬ 
able. As seen see Figs. ,b, a macroscopic fraction of 
the population synchronously follows persistent orbital 
trajectories with a constant radius, see also Supplemen¬ 
tary Movie 2 [IB] . 

This alternative type of collective motion is well cap¬ 
tured by two order parameters: the global polarization 
n and the global spin S, which is defined as the absolute 
value of the most probable instantaneous average spin 
{o’i{t))i. These two quantities sharply increase as rj goes 
below as shown in Fig. [^. It is also worth noting 
that, the temporal variations of (o’i(t))i is intermittent 
over time scales that are much larger than the individ¬ 
ual relaxation time t regardless of the noise amplitude, 
see Fig. The rare intermittent events correspond to 
the interruption of the circular trajectories. The parti¬ 
cles randomize their orientation at the same time and 
subsequently self-organize again to rotate synchronously, 
possibly along a different direction. This long-time dy¬ 
namics results in a bimodal structure of 'P{{ai{t))i), seen 
in Fig. and a posteriori justifies defining the spin or¬ 
der parameter as the modulus of the most probable value 
of E(f). 

In addition the individual circular trajectories ob¬ 
served in this spinning phase are not robust: they con¬ 
tinuously form and disrupt, and neither H nor S con- 



FIG. 2: Color online. Case a = ^. (a)Typical trajectories in 
the spinning phase ( 7 ; = 4 x 10”“^). (b) Snapshot of the posi¬ 
tion and orientation of the motile particles for the same noise 
value. The color code accounts for the instantaneous spin as 
in Fig. 1 (j; = 4 X 10“"^). (c) Variations of the order pa¬ 
rameters n (light/black symbols) and E (dark/red symbols) 
with respect to the noise amplitude rj. (d) Time series of the 
instantaneous global spin E(t), and (e) related probability 
distribution that displays a clear bimodal structure, (f) The 
spatial condensation at the transition is clearly seen from the 
variations of Sr with rj, the horizontal line corresponds to the 
value of Sr for a uniform distribution of particles in a square 
box. All the simulations correspond to p = 8, uo = 0.5 and 
L = 16 (A = 2048). 


verges to 1 as ry goes to 0 (Fig. [!=)■ At any given time, a 
small yet finite fraction of particles does not follow syn¬ 
chronous circular trajectories even for vanishingly small 
noise. These unusual strong fluctuations in the weak- 
noise limit, together with the circular motion of the par¬ 
ticle, will be better understood in the next section, when 
inspecting the phase diagram of the general a model. 

Last, we stress that at the onset of global spinning 
motion, the spatial distribution of the particles is hetero¬ 
geneous. The population concentrates into a large denser 
region surrounded by a more dilute ensemble of spinning 
particles. This spontaneous breaking of the translational 
invariance is quantified by measuring the variations of 
Sr = |(ri — (ri)i)i_t| as a function of the angular noise. 
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FIG. 3: Color online. Alignment and anticipation for finite a. (a) Phase diagram showing the region where the isotropic, the 
flocking and the spinning phases are stable, (b) Sharp variations of the polarization and spin order parameters with rj (same 
color code as in Fig. [^. (c) Exemples of particle trajectories for a = 0 . 87 r /2 in the flocking state, for rj = 10“^. (d) Variations 
of the angle ip between the polarization and the principal axis of the simulation box, as a function of the anticipation angle 
OL. (e) Typical spinning trajectories observed for 77 = 7 x 10“^. (f) Instantaneous position and orientation of one half of the 
particles, with the same color code for the spin variable. The entire population rotates in the same direction, (g) Variation of 
the curvature of the spinning trajectories with the noise amplitude for a = 0 . 47 r. (h) Variation of the curvature of the spinning 
trajectories as a function of a in the zero-noise limit. Full line: analytical prediction k = sina/uor. For all simulations, p = 8 , 
Vo = 0.5, and L = 16 (iV = 2048). 


Fig. [2f . This spatial heterogeneity emerges in the ab¬ 
sence of explicit couplings between the positional degrees 
of freedom, which is a very persistent feature of all the 
active-particle models involving any form of short-range 
interactions [19]. 

C. Case a > 0: From flocking to spinning and 
swarming. 

We now discuss the general case, for which the antic¬ 
ipation angle a can take any positive value. The col¬ 
lective behavior of the population is summarized by a 
phase diagram that divides the ( 77 , a) plane into three 
regions as shown in Fig. The boundaries of the phase 
diagram correspond to the points where the order pa¬ 
rameters, n or E, go from zero to a finite value. When 
0 < a ^ 7 r/ 2 , for very high noise values the system forms 
an isotropic and homogeneous gas. Reducing 77 , the pop¬ 
ulation undergoes two subsequent phase transitions to¬ 
ward collective motion. A flocking transition occurs at 
V = tigf{o)'- the polarization increases sharply while the 
global spin remains vanishingly small. Fig. However, 
unlike all the Vicsek-like models with periodic boundary 


conditions, the global polarization does not align with 
one of the principal axis of the simulation box. In con¬ 
trast, it makes a finite angle, ip, that is set by the antic¬ 
ipation angle a, Figs. [^, and p. 

When further reducing the noise amplitude below 77 = 
VFsicK), the flock self-organizes into a spinning phase. As 
opposed to the case where a = 7r/2 examined in previous 
section, we show in Fig. P that E increases sharply and 
saturates to 1 together with H, as 77 goes to 0 deep in 
the spinning phase. The individual circular trajectories 
are now robust. They do not disrupt and reform inter¬ 
mittently; all the particles endlessly follow synchronous 
circular orbits, see Figs. |^, and [^. 

In the spinning phase, the curvature k of the particle 
orbits has a constant value, Fig.j^. This result can be ex¬ 
plained by considering the zero-noise limit. When rj — 0, 
a perfectly polarized state with all particles aligned is a 
solution of Eq. for any value of a, provided that the 
particles rotate at the same angular velocity 9 = sin a. 
As the particle velocity Vq is a constant, the correspond¬ 
ing trajectories are circles of curvature k = sina/rioT. 
As seen in Fig. P, this mean-held argument correctly 
accounts for our numerical results when 0 < a < 7 r/ 2 . 

Figure P shows that the phase behavior of the pop- 
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ulation is qualitatively different when a approaches and 
exceeds 7 r/ 2 . The polar-liquid/flocking phase, does not 
exist anymore above ax 0 . 9 ( 7 r/ 2 ). The system under¬ 
goes a single transition from an isotropic to a spinning 
state. For ax ~ 0 . 9 ( 7 r/ 2 ) the transition occurs through a 
tricritical point, where the flocking, the spinning and the 
isotropic phases coexist. In addition, above ax the spin¬ 
ning phase displays marked qualitative differences from 
the one characterized above for small values of a: the 
polarization and the spin do not saturate to 1 as ry —>■ 0 
(as previously noted for a = 7 r/ 2 ). In addition the curva¬ 
ture of the trajectories deviates from the naive mean-field 
prediction. Fig. ii- 

In order to elucidate this rich phase behavior, we now 
investigate the linear stability of the two ordered states. 
Let us first consider an ensemble of particles forming a 
polar liquid oriented along 0 = 0. Far from a transition 
point, the angles weakly deviate from 0 ^ = 0 and the 
angular dynamics reduces to: 

0i(t) -- sina((Ti)jga, - cosa(0i - 0j)jGni + (4) 

In this polarized but non-spinning state, the individual 
spins undergo uncorrelated fluctuations. Therefore, for 
any finite interaction radius, the first term on the right- 
hand side acts as an additional angular noise that im¬ 
pedes the velocity alignment, in stark contrast with the 
intuitive picture that one could have about the effect 
of anticipation. The magnitude of this effective angu¬ 
lar noise increases with a which qualitatively explains 
why the flocking transition occurs for smaller values of 77 
as a increases. In addition, it readily follows from Eq. 
that the effective angular stiffness, cos a, is negative for 
a > 7 r/ 2 . Hence the homogeneous polar-liquid phase is 
linearly unstable to angular fluctuations. In agreement 
with the phase diagram shown in Fig. 3a, no polar-liquid 
phase exists for a > 7r/2. 

Let us now repeat the same analysis for the spin¬ 
ning phase. Introducing the angular fluctuations 50^ = 
Oi zLt sin a deep in the ordered phase, the orientational 
dynamics reduces to: 

69i{t) - cos a{S8t - S8j)jen, + (5) 

Again we find that the orientational fluctuations are sta¬ 
bilized for a < 7 r /2 and amplified otherwise. This stabil¬ 
ity analysis is consistent with the strong fluctuations of 
the individual trajectories found for a = 'kI2\ the spin¬ 
ning state is linearly unstable for large anticipation an¬ 
gles. The ordered spinning phase numerically found for 
a > 7r/2 is therefore stabilized by the nonlinear nature 
of the spin-angle interactions and cannot be captured by 
a linear-stability analysis alone. 

The flocking and the spinning phases are separated 
by another dynamical state where the polar-liquid phase 
condenses into compact swarms undergoing coherent di¬ 
rected motion. This behavior is characterized and illus¬ 
trated in Figs. 1^ and respectively, see also Sup¬ 
plementary Movie 4 m- This swarming phase exists 



FIG. 4: Color online. Swarming at the onset of collective 
spinning, (a) The variations of Sr echo the condensation of 
the particles into a single compact swarm undergoing directed 
motion, (b) Snapshot of the instantaneous particle positions 
and orientations in the swarming state where If = 0.977 and 
E = 0.063. a = 0 . 877 / 2 , T/ = 3 X 10“^, p = 8 , uo = 0.5 and 
L = 16 (A = 2048). 


in a narrow range of the noise amplitude 77, typically a 
few percent above r/ps- Unlike the band-shape patterns 
observed in the standard Vicsek models, these compact 
swarms are not surrounded by a sea of isotropic parti¬ 
cles, but freely propagate in an empty space. To our 
knowledge, this condensation is the first evidence of the 
self-organisation of motile particles into compact swarms 
without any attractive interactions. 

We close this section by pointing two questions that 
readily arise from the above analysis. First the nature 
of the transition between the four phases (gas, flocking, 
spinning and swarming states) remains to be elucidated 
(critical, first order, or sharp cross-over). Answering this 
question would require studying much larger systems and 
a systematic finite-size analysis which goes beyond 
the scope of the present paper. A second natural question 
concerns the robustness of these phases with respect to 
the noise structure. We focused here on the case of an 
angular noise. Whereas the present results are expected 
to hold for a vectorial noise as well, the impact of more 
complex multiplicative noises, such as the one found to 
be relevant to locust swarms in |22j . remains an open 
question. 


IV. CONCLUSION AND PERSPECTIVES 

We conclude this paper by addressing potential ap¬ 
plications of our hndings to multi-agent robotics. In 
this context, biomimetic strategies have always been ap¬ 
pealing to achieve collective intelligence. However man¬ 
made programable units could exploit a virtually infinite 
number of alternative strategies to accomplish emergent 
tasks. In this paper, we have demonstrated that sim¬ 
ple artificial interaction can be effectively used to self- 
organize an entire population. By coupling the spin and 
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the orientation of motile particles, three symmetries can 
be spontaneously broken (rotation, translation and spin- 
reversal) thereby resulting in three spatiotemporal pat¬ 
terns: homogeneous and coherent directed motion (flock¬ 
ing) , synchronous circular motion (spinning) and conden¬ 
sation into a compact group propagating in a coherent 
fashion (swarming). In addition, the transition from a 
swarming to a spinning state can be triggered by a minute 
increase of the anticipation angle a that characterizes the 
spin-orientation coupling. Therefore this sharp transition 
could provide a useful design rule to arrest and disperse 
a motile swarm without having to stop the individual 
propulsion. Since all the particles undergo synchronous 
circular orbits in the spinning state, a swarm of particles 


moving coherently could stop and explore the neighbor¬ 
hood of a given spot without having to worry about pos¬ 
sible collisions between the motile units. Experimental 
investigations along these lines are now needed to ad¬ 
dress the relevance of these simple interaction rules to 
devise functional robotic swarms. 
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